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We study pattern formation in the bounded confidence model of opinion dynamics. In this 
random process, opinion is quantified by a single variable. Two agents may interact and reach a 
fair compromise, but only if their difference of opinion falls below a fixed threshold. Starting from 
a uniform distribution of opinions with compact support, a traveling wave forms and it propagates 
from the domain boundary into the unstable uniform state. Consequently, the system reaches a 
steady state with isolated clusters that are separated by distance larger than the interaction range. 
These clusters form a quasi-periodic pattern where the sizes of the clusters and the separations 
between them are nearly constant. We obtain analytically the average separation between clusters 
L. Interestingly, there are also very small quasi-periodic modulations in the size of the clusters. The 
spatial periods of these modulations are a series of integers that follow from the continued fraction 
representation of the irrational average separation L. 

PACS numbers: 89.75.Kd, 82.40.Ck, 05.45.-a 


The so-called “bounded-confidence” model [IH3 and 
variants thereof have been widely used to model opin¬ 
ion dynamics dSl and have attracted a considerable 
amount of interest . The bounded confidence model is 
appealing because it captures the tendency for reaching 
compromise through social interactions, while also taking 
into account a certain degree of conviction. Numerical 
studies show that political parties emerge in the bounded 
confidence model as a result of a pattern formation pro¬ 
cess [3, Si- In this letter, we obtain analytically the 
wavelength governing this process. 


We focus on a version of the bounded confidence model 
in which opinions are quantified as discrete variables 
1 < n < A^. In each interaction, two agents with opin¬ 
ions ni and n 2 change their initial opinions by adopt¬ 
ing the average opinion (ni,n 2 ) ^ such 

a compromise occurs only when the opinion difference is 
smaller than some fixed threshold |ni — n 2 | < cr. We 
set the threshold a = 2 and exclude interactions between 
agents whose opinion difference equals one, \ni —n 2 \ = 1, 
to ensure that opinions remain discrete variables. In this 
simplified version of the bounded confidence model, opin¬ 
ions change according to Q 


(n — 1, n -I- 1) —)• (n, n). 


( 1 ) 


Clearly, this process conserves population and opinion. 

Let Pn{t) be the probability density of agents with 
opinion n at time t. This density obeys the rate equation 

dP 

, = ‘2‘Pn-lPn+l — Pn{Pn-2 + Tri-|-2)- (2) 

dt 

In writing this equation, we implicitly take the infinite 
population limit. It is simple to check that Eq. m con¬ 
serves population, and opinion, '^^nPn- 

The initial distribution of opinions is uniform with 


compact support, 

{ 0 u < 1, 

1 l<n<A, (3) 

0 N < n. 

We view the parameter N as the opinion “spectrum,” 
and also note that N is the only parameter in the model. 
The evolution equation ([2]) is invariant under the scaling 
transformation P —>■ aP and t —> t/a and hence, we may 
set the uniform initial density to unity. This choice allows 
us to compare systems with different opinion spectrums. 

The nature of the interaction o, also reflected by the 
evolution equation ([2]) , implies that any probability den¬ 
sity that satisfies Pn-iPn+i = 0 for all n is stationary. 
Clearly, there are infinitely many such steady-state solu¬ 
tions. Starting from the (unstable) initial condition ([3]), 
the deterministic rate equation (ED evolves the system 
toward one of those (stable) steady-state solutions 0. 

In the final state (i.e., in the limit t —>■ oo), the sys¬ 
tem reaches a steady state where P„_i(oo)Pti+i(oo) = 0 
for all n. In this state, there are multiple opinion “clus¬ 
ters” with each cluster localized to two neighboring lat¬ 
tice sites (see figure [1]). These clusters are noninteracting 
because the separation between them exceeds the inter¬ 
action range. 

To quantify the size and opinion of each cluster, we 
compute for every pair of occupied lattice sites the 
mass m = P„(oo) -I- P„+i(oo) and the non-integer posi¬ 
tion X = [nP„(oo) -I- (n -I- l)P„+i(oo)]/m. Let be the 
mass of the ith cluster and Xi be the position of the ith 
cluster. Conservation of population and opinion sets the 
sum rules mi = N and J2i = N{N l)/2. 

Unlike the probability density, the cluster mass forms a 
quasi-periodic pattern (see figure [T]) as clusters of nearly- 
identical mass are separated by a nearly-identical dis¬ 
tance. This pattern can be characterized by the average 
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FIG. 1: [bottom] The probability density P„(oo) versus n. 
[top] The cluster mass m versus position x. . Shown are 
results of integration of ©-([S]) with N = 109. The dashed 
line corresponds to the theoretical value L — 5.671820 


separation L between clusters. 


L = lim (to). 

W-s-oo 


(4) 


With this definition, the average number of clusters scales 
as N/L in the limit iV —>■ oo. Previous numerical stud¬ 
ies reported the value L 5.67 0. In this letter, we 
use theoretical methods to analyze the evolution of the 
probability density Pn and analytically obtain L as the 
wavelength that governs the underlying pattern forma¬ 
tion process. 

The uniform initial state (jS]) is unstable with respect to 
perturbations that propagate from either boundary into 
the unstable uniform state 15, in. By substituting the 
small periodic disturbance 


Pn{t) — 1 oc exp[i(A:n — ojt)] (5) 

into the evolution equation © we find the dispersion 
relation between frequency uj and wavelength k, 


w = 21(2 cos/c — cos 2fc — 1). (6) 


Because the quantity —iuj is positive for 0 < fc < 7r/2, 
perturbations with wavenumber in that range initially 
grow exponentially with time. The fastest growing mode, 
by ordinary linear stability analysis, follows immedi¬ 
ately from ([6]). The maximum of —iuj in ([6]) is set by 
dhj/dk = 0 which yields linear = 7r/3 or alternatively 

Tlinear — ^Tt/ /^linearj that IS, Tiinear — 6. 

The perturbations propagate from the stable state into 
the unstable state at a constant velocity v (see figure 
[21). A saddle point analysis shows that the propagation 
velocity v obeys (for a comprehensive review see [l7j| and 
also [isl - Elj ]! 


dhj Im[w] 
dk Im[fc] 



FIG. 2: Propagation from the stable into the unstable state. 
Shown is Pn{t) at equally spaced time intervals versus n. The 
curves Pn{t) are shifted vertically with the earliest time at 
the bottom and the latest time on top. 


The solution to this equation is the complex wavenumber 
/c* = /cfront + with /cfront = 1.183032 The constant 
A = 0.467227 characterizes the exponential decay of these 
periodic perturbations — 1 ^ 

The wavelength of perturbations at the leading edge of 
the front is Tfront = 27r/fcfront or explicitly, 

Tfront = 5.311086. (8) 

Moreover, the propagation velocity v = Im[r(;,]/Im[A:*] 
where w* =w{k^) is 


V = 3.807397. 


(9) 


Our numerical results confirm that in the leading edge 
of the propagating front, the wavelength of the periodic 
deviations from the uniform state is indeed given by ([8]) . 

Far behind the traveling wave, that is, in the wake 
of the wave, the system reaches a steady state with 
Pn-iPn+i = 0 for all n <^vt. In this region, clusters are 
fully-developed. Interestingly, far behind the propagat¬ 
ing front, the pattern that forms, and which ultimately 
controls the spacing between clusters, has a larger wave¬ 
length due to a Doppler-like effect. The frequency of 
oscillations in the front w* — k^,v, measured in the co¬ 
moving frame, translates to a zero frequency in the rest 
frame and hence, to stationary patterns, precisely for the 
wavenumber TZ, 


= ( 10 ) 

V 

We note that unlike the Doppler effect, the normalized 
shift in wavenumber (fc* — does not equal a ratio of 

two velocities as it is complex. The resulting wavenumber 
m is k = 1.107789 and the corresponding wavelength 
L = 2Tr/k is 


V = 


(7) 


L = 5.6718200283. 


( 11 ) 
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FIG. 3: The cluster mass m versus position x. The red line 
corresponds to the theoretical value dl) and the blue line, to 
a 100-point running average. 



FIG. 4: The front location Xf versus time t. 


Hence, out of the entire range of possible wavelengths cor¬ 
responding to linearly unstable perturbations, 0 < L < 6, 
the wavenumber m is “selected” by the dynamics of 
Eq. Q. We also note the inequality Lfront < L < Linear ■ 
Our numerical results give excellent confirmation of the 
theoretical prediction (see Fig. El): the numerically- 
measured wavelength L = 5.67185 is within 10“® of (ED- 

To efficiently perform the computation, we integrated 
the equations using a lattice with fixed size J\f that is 
moving at the same speed as the traveling wave. Nu¬ 
merical integration in this co-moving reference frame is 
feasible because far ahead of the traveling wave = 1 
and far behind it, the system settles into a steady-state. 
The lattice is shifted by one lattice site, n ^ n+1 when¬ 
ever the deviation from the uniform state at the extreme 
lattice site, far ahead of the traveling front, exceeds an 
infinitesimal threshold, \Pjg- — 1| > e. We used a Runge- 
Kutta (4,5) integration method with adaptive step size 
below 10“^ on a lattice of size Af = 2, 000 and the thresh¬ 
old e = 10“^®°, well above the smallest machine precision 
10“308^ To resolve the leading edge at this precision, we 
integrated the equations for = P„ — I where the lead¬ 
ing edge decays to the constant state Qn = 0. Whenever 
the lattice is shifted by one site, time and the probabil¬ 
ity density at the site n = 200, well behind the leading 
edge, were recorded. This approach allows us to integrate 
the equation to times t « 3 x 10® and effectively, study 
systems with very large opinion spectrums N si 10®. 

The number of shifts at time t directly measures the 
front position Xf. As shown in figured The numerically- 
measured propagation velocity v = 3.80732 is in excel¬ 
lent agreement with the theoretical prediction It 

is remarkable that our computation, for which the ra¬ 
tio between system size and wavelength is moderate, 
Af/L « 350, yields such high-precision measurements of 
L and v. In general, a cutoff error e in the propagat¬ 
ing front results in logarithmic correction S (lne)“^ in 
the propagation velocity [1^. Moreover, the cutoff de¬ 
cays exponentially with system size, e ^ exp(—AA/”). By 
combining these two scaling laws, we find the algebraic 


relationship between system size Af and correction 6 to 
the velocity, <5 ^ {XAf)~^. The velocity correction we 
observe, S ~ 10“® for Af ~ 10®, is consistent with this 
scaling law. 

The irrational wavelength L in m is not commensu¬ 
rate with the unit lattice spacing. As a result, the pat¬ 
terns are not strictly periodic but rather, they are quasi- 
periodic (see hgureE])- Interestingly, we observed small 
but striking “super-patterns” induced by near-resonances 
between L and the lattice spacing. The integer periods of 
these super-patterns are found from a continuous-fraction 
expansion of the wavelength 


L = 5 


1 + 


= 6 . 


17 363 1469 


2 - 1-7 


3 ’ 64 ’ 259 


( 12 ) 


Indeed, 3 clusters can be accommodated within the in¬ 
teger period 17. Figure [2] for the probability density P„ 
shows that a large peak in the quantity is usually fol¬ 
lowed by two smaller ones. Hence, the system exhibits 
quasi-periodic behavior with integer period 17. 

Furthermore, there is also quasi-periodic arrange¬ 
ment of clusters with integer period 363. Figure 
El shows that the cluster mass varies in the range 
L — A < m < L + A. The variation in cluster mass is 
very small A/m ~ 2 x 10“®. As a function of position, 
these small variations in cluster mass repeat with integer 
period 363. According to the continued fraction (I12|) . 64 
clusters form this pattern. This intriguing behavior was 
overlooked in previous studies that used much smaller 
values of the opinion spectrum A [J [^ . 

Interestingly, the simple evolution equation ([2]) leads to 
a hierarchical patterns, governed by a series of integer pe¬ 
riods. The primary pattern, as shown in figure [1] consists 
of a nearly-periodic arrangement of clusters with nearly- 
identical separation. On the first hierarchical level. Mi 
clusters are arranged in a nearly periodic super-pattern 
with period Li. On the second hierarchical level, M 2 
patterns form a more intricate super-pattern with the 
integer period L 2 - The fractions Mi/Li, M 2 /L 2 , and so 











FIG. 5: The cluster mass m versus position x (large square) 
and versus shifted position x + 363 (small circle). The line 
corresponds to the theoretical value m- 
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FIG. 6: The cluster mass m versus position x (large square) 
and versus shifted position x + 1445 (small circle) for the 
three-site interaction model in dSJ. The line corresponds to 
the theoretical value m- 


on are rational approximations of the wavelength L that 
follow from its continued fraction representation (HU. 

To examine the robustness of the above behavior, we 
also considered the compromise process (n — 1, n -|- 2) — 
(n,n-\-l) where the interaction range is three lattice sites. 
In this case, the probability density evolves according to 
the rate equation 

dP 

= Pn-2Pn+l + Pn-lPn+2 ~ Pn{Pn-Z + ^rt+3)-(13) 

This equation conserves population and opinion. Start¬ 
ing from the initial condition ([3]) , the probability density 
evolves toward a steady state where Pn(oo)Pn+ 3 (po) = 0 
for all n, and therefore, opinion clusters are now localized 
to three consecutive lattice sites. By substituting ([5]) into 
(d, the dispersion relation is 

w = 2i(cosfc-I-cos2fc — cosSfc — 1). (14) 

By repeating the analysis leading to (d) we obtain the 
average separation between clusters and propagation ve¬ 
locity 

L = 8.5502770500 and r; = 2.50631. (15) 


Also, the wavelength of patterns nucleating at the front is 
Afront = 8.02282. Numerical integration of the evolution 
equation gives L = 8.5503 and v = 2.5063, in excellent 
agreement with the theoretical predictions. 

Figure |6] demonstrates the emergence of super¬ 
patterns. In this case, continued fraction representation 
of the wavelength d gives the rational approximations 


L = 


1 -h 


= 9 


17 77 171 1445 


1 +^ 


2 ’ 9 ’ 20 ’ 169 


.(16) 


'2TW 


Figure |6] shows that modulations in cluster mass are pe¬ 
riodic and well-characterized by the integer period 1445, 
that follows from the continued fraction d- Accord¬ 
ingly, the super-pattern consists of 169 clusters. Fur¬ 
thermore, the amplitude of the variations is very small. 
A/m ss 5 X 10-4. 


We also examined a linear interpolation between (P|) 
and d, modeling a compromise process where second¬ 
er third-neighbor interactions occur with relative weights 
r and 1 — r, respectively. We found excellent agreement 
between the theoretical predictions and the numerical re¬ 
sults for the wavelength L and the velocity v for all values 
of T. As expected, the wavelength decreases monotoni- 
cally with the mixing parameter r, which decreases the 
effective interaction range. Surprisingly, however, the ve¬ 
locity is not a monotonic function as it reaches a mini¬ 
mum for r Ri 0.88. Thus, replacing a small fraction of 
the second with third-neighbor interactions slows down 
the spreading of the compromise process. 

Moreover, we do not observe “mode-locking”. In this 
phenomenon, common in pattern forming systems that 
are exposed to an external spatially-periodic forcing, 
there is pronounced locking of the wavelength to that 
of the forcing 2J, [ 23 . In lattice systems, one there¬ 
fore expects the observed wavelength to be an integer, 
in resonance with the unit lattice spacing, whenever the 
predicted wavelength is close to an integer. In our case, 
the predicted wavelength L is integer (6 or 7 or 8) for 
particular values of the mixing parameter r. We studied 
values of r near those resonances but did not observe any 
locking: the wavelength varies smoothly and adheres to 
the predicted values. Further, we did not observe sub¬ 
harmonic resonances as in [26l429l| . where the observed 
wavelength is an integer multiple of the predicted wave¬ 
length. 

We now briefly discuss the original bounded confidence 
model introduced by Weisbuch et al [ij. In that model, 
opinion is quantified by a continuous variable 0 < x < N 
with N the opinion spectrum. Agents can interact and 
reach fair compromise but only if their opinion difference 
falls below a fixed threshold, set to unity without loss of 
generality Q 

if \X1-X2\<1. (17) 

The probability density P{x,t) of agents with opinion x 
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at time t obeys the evolution equation Q 


-P{x,t) = jj dxidX 2 P{xi,t)P{X 2 -,t) 


|a^l-a:2|<l 


5 { X — 


Xi + X2 


— 6{x — xi) 


(18) 


If the restriction of the integration range is ignored, this 
equation describes inelastic collisions [30|-l33|. According 
to the interaction (13, opinion clusters are now perfectly 
localized (delta-functions) and in the steady-state these 
localized clusters are separated by distance larger than 
unity. 

Consider the uniform initial condition: P{x, 0) = 0 
for a; < 0 or a; > and P(a:, 0) = 1 for 0 < x < iV. 
This state is unstable with respect to perturbations that 
propagate from the boundary into the unstable uniform 
state. According to (fT51) . a small periodic disturbance 
P{x, t) — 1 (X exp[*(fcx — wt)] has the dispersion relation 


uj = 2i 


O sin(fc/2) 
k/2 



(19) 


The fastest growing mode follows from dui/dk = 0 which 
yields /cunear = 2.7906 or alternatively Tunear = 2.2515. 
The solution to © is now fc* = fcfront + *A with 
fcfront = 3.083750. The decay constant A = 1.294620 
characterizes the exponential decay far into the un¬ 
stable state, (j){x) ^ exp[—A(x — vt)]. The wave¬ 
length of perturbations at the propagating front is 
Tfront = 2.037514. The propagation velocity v = 
Im[w*]/Im[A;*] is u = 0.794754. Far behind the propagat¬ 
ing front, that is x <C vt, localized clusters form, and 
these clusters are separated by distance L. The cor¬ 
responding wavenumber is fc = 2.924255 and the corre¬ 
sponding wavelength is 


L = 2.1486444707. (20) 


The wavelength estimated using numerical integration re¬ 
sults for relatively small values of Q, L ss 2.155, is 
reasonably close to the exact result (1^ . 

In summary, we studied pattern formation in the 
bounded confidence model of opinion dynamics. Our fo¬ 
cus was the wavelength that governs the mosaic of frozen 
clusters that develop, starting from a uniform state. We 
obtained analytically the two wavelengths that govern 
the pattern formation process: the wavelength of pertur¬ 
bations at the leading edge of the traveling wave front and 
the wavelength of the resulting patterns in the wake of 
the wave. We examined discrete and continuous versions 
of the bounded confidence model. In the former case, we 
verified the theoretical predictions using high-precision 
numerical measurements of the pattern wavelength and 
propagation velocity. 

The wavelength of the patterns is irrational and since 
it is not commensurate with the regular lattice, the pat¬ 
tern formation process is hierarchical. Frozen clusters 


constitute the “building bocks” in this hierarchy. Integer 
number of clusters form quasi-periodic structures and the 
period of these super-patterns is an integer, too. Next, 
a larger number of clusters form a more intricate super¬ 
pattern with a larger integer period. The numbers of 
clusters and the periods that characterize these super¬ 
patterns follow from continued fraction representation of 
the irrational wavelength governing the pattern forma¬ 
tion process. 

We observed that not all rational approximations of 
the wavelength necessarily correspond to a super-pattern 
(see figure [6|). Further analysis is therefore needed to un¬ 
derstand why certain integer fractions are realized while 
others are not, and more generally, to characterize the 
intricate structures of the super-patterns. 
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